Singular limit of a two-phase flow problem in porous 
medium as the air viscosity tends to zero* 



R. Eymard ^, M. Henry ^ and D. Hilhorst ^ 
Abstract 

In this paper we consider a two-phase flow problem in porous media and study 
its singular limit as the viscosity of the air tends to zero; more precisely, we prove 
the convergence of subsequences to solutions of a generalized Richards model. 

1 Introduction 

Hydrologists have studied air-water flow in soils, mainly using the so-called Richards 
approximation. At least two hypotheses are physically required for this model to be 
applicable: the water pressure in the saturated region must be larger than the atmospheric 
pressure and all the unsaturated regions must have a boundary connected to the surface. 
However, in many situations, these hypotheses are not satisfled and a more general two- 
phase flow model must be considered. This work explores the limit of this general model 
as the viscosity of the air tends to zero, which is one of the hypotheses required in the 
Richards model. To that purpose we prove the existence of a weak solution of the two- 
phase flow problem and prove estimates which are uniform in the air viscosity. In this 
paper, we assume that the air and water phases are incompressible and immiscible. The 
geometric domain is supposed to be horizontal, homogeneous and isotropic. Our starting 
point is the following two-phase flow model, which one can deduce from Darcy's law 

,^ . j ut- div{h^{u)V {p)) = 

^ {l-u)t-dzv{j^ka{u)V{p + Pc{u))) = Sa, 
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where u and p are respectively the saturation and the pressure of the water phase, kyj 
and ka arc respectively the relative permeabilities of the water and the air phase, /i is 
the ratio between the viscosity of the air phase and that of the water phase, is the 
capillary pressure, s^, is an internal source term for the water phase and Sa is an internal 
source term for the air phase; these source terms are used to represent exchanges with the 
outside. We suppose in particular that the physical functions k^, ka and Pc only depend 
on the saturation u of the water phase, and that kw{l) = ka{0) = 1. The aim of this 
paper is the study of the limit of the two-phase flow problem as // J, 0. 



The classical Richards model as formulated by the engineers is given by 

(7^) 



U=p-\patm -P)- 



where the properties of capillary pressure Pc — Pc{u) are describes in hypothesis (Hs) 
below. For the existence and uniqueness of the solution of Richards model together with 
suitable initial and boundary conditions as well as qualitative properties of the solution 
and methods for numerical approximations we refer to [1], [6], [10], [11]. In this article, 
we will show that the singular hmit as | of the two phase flow problem (TV) has the 
form 

Ut — div{ku,{u)Vp) = Sw 

M = 1 or V{p + Pc{u)) = a.e. in x (0,T). 
We remark that a solution of (Tl) with u> satisfies [TBV) . 



{TBV) 



This paper is organized as follows. In Section 2 we present a complete mathematical for- 
mulation of the problem, and state the main mathematical results, which include a precise 
formulation of the singular limit problem. We give a sequence of regularized problems 
in Section 3, and prove the existence of a classical solution. In Section 4 we present a 
priori estimates, which are uniform in an extra regularization parameter 5 and in the air 
viscosity /i. In Section 5, we let 5 J, and prove that the solution converges to a solution 
of the two phase flow problem. We study its limiting behavior as the air viscosity // tends 
to zero in Section 6. Finally in Section 7 we propose a finite volume algorithm in a one 
dimensional context and present a variety of numerical solutions. 
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2 Mathematical formulation and main results 



We consider the two-phase flow problem 

(1 - u)t = div(^-ka{u)V{p + Pc{u))j 



+ {1 - f'^{c))s - {1 - nu))s, 

p{x, t)dx = 0, 

Vp.n = 0, 
V{p + Pc{u)).n = 0, 
, u{x,0) = uo{x), 

where T is a positive constant, Qt :— fl x (0, T) and where we suppose that 



in Qt, 


(2.1) 


in Qt, 


(2.2) 


for t e (0,T), 


(2.3) 


on dn X (0,T), 


(2.4) 


on on X (0,T), 


(2.5) 


for X E il. 


(2.6) 



{He 



Q is a smooth bounded domain of where the space dimension N is arbitrary, 

Um e (0, 1), 

ceL°°{Qx (0, T)) and u^<c<l, 
uq G L°°{Q) and Um < Uo < 1, 

s e L2(fi), s > 0, se L'^iVt), s > and J {s{x) - s{x))dx = 0, 

G C\[0, 1]), k'^ > 0, fc^(O) = 0, kM = 1 and k^uj > 0, 
fca G C2([0, 1]), A;;<0, A;„(1)=0, A;„(0) = 1 and fc„(s) > for all s G [0, 1), 
Pe e C0([0, 1]) U C^iO, 1)), < and sup,,[o,i)(-^a(«K(s)) < +oo, 
/xG(0, 1]. 



In this model, u and p are respectively the saturation and the pressure of the water 
phase, kyj and ka are respectively the mobilities of the water phase and the mobility of 

the non-water phase and p,. is the capillary pressure. We assume in particular that the 
permeability functions fe^,, ka and the capillary pressure pc only depend on the saturation 
u of the water phase. Here, we suppose that the flow of the water phase in the reservoir 
is driven by an injection term /^(c)s and an extraction term f^{u)s where s and s are 
given space dependent functions, c is the saturation of the injected fluid; if c = 1, only 
water will be injected, if c = 0, only air will be injected, whereas a mixture of water and 
air will be injected if < c < 1. The function is the fractional flow of the water phase, 
namely 



kw{s) 
M'^(s) 



, with M^(s) 



kw{s) + -ka{s). 
A* 



(2.7) 



In particular, we remark that 



f^{s) is non decreasing. 



(2i 
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Next we introduce a set of notations, which will be useful in the sequel. 

9{S) ^ - r ka{T)p',{T)dT, (2.9) 

Jo 



as) = ^ka{T)p',{T)dr, (2.10) 

Q'is)= r nr)Pc{r)dr, (2.11) 
Jo 

and 

for all s e [0, 1] . This implies in particular that 

n>'{s) + Q'^(s) = pc{s) - pe(0), for all s e [0, 1]. (2.13) 
Definition 2.1 The pair {u^,p^) is a weak solution of Problem (S^) if 
G L°°((] X (0, T)), with < m'^ < 1 'm Qt, 

p'^ e L2(0, T; //HJl)), / p^lx, t)dx = /or almost every t e (0, T), 

^(ii'^)eL2(0,r;i/i(Q)), 

r / u^'iftdxdt^ r [ k^{u^')Vp^.Vipdxdt- r [ (f^{c)s- f^'{u^)s\ipdxdt 
Jo Ja Jo Jn Jo Jn \ / 

- / Uo(x)(p(x,Q)dx, (2.14) 
Jn 

and 

u^^iptdxdt ^ £ ^kaiu^") (vp" + Vpc(m'')) S/ipdxdt 
- |J_£((l-r(c))s-(l-r(M'^))s)<^dxdt-^(l-«o(a;))9:'(a;,0)cia;, (2.15) 
for all ip inC:^{w e W2'\Qt),w{., T) = m Q}. 

Our first result, which we prove in Section 3, is the following 

Theorem 2.2 Suppose that the hypotheses (ifi) — (Hg) are satisfied, then there exists a 
weak solution {u'^,p^) of Problem (5*^). 



Next we define the discontinuous function x by 

' if s e [0, 1) 



as well as the graph 



His) := 



1 if s = 1, 

if s e [0, 1) 



[0, 1] if s = 1. 

The main goal of this paper is to prove the following convergence result, 

Theorem 2.3 Suppose that the hypotheses {Hi) — {Hq) are satisfied, then there exists a 
subsequence {{u^" , p^"))neN of weak solutions of Problem {S^") and functions u, p, f such 
that 

u e L^{Qt): <u<1 rnQx, 
feL^iQr), 0</<l inQr, 
peL'{0,T;H\n)), 
ka{u)Vpc{u) e L-\n X {0,T)), 

and 

{u'''")neN tends to u strongly in L^{Qt), 
{p'^")neN tends to p weakly in L'^{0,T; H^{fl)), 

as /In tends to zero and 



JQ 



uiptdxdt 



/ kyj{u)'Vp.'V(pdxdt — / / [x{c)s — fs](pdxdt 
in Jo Jn\ J 



— / uo(x)(p(x, 0)dx, 
Jn 



(2.16) 



for all (fi & C, where f{x,t) G H{u{x,t)) for {x,t) e Qt- Moreover we also have that 



T 

JQ, 



ka{u) 



1 2 



Vp + Vpc{u) 



dxdt — 



and 



/ p{x,t)dx = 0, for almost every t G (0,T). 
Jn 

Formally, u satisfies the following limit problem 

ut = div[ky,{u)Vp^ + x(c)s - /s, in Qt, 

' Vu.n = 0, on x (0,T), 

, u{x, 0) = Uq{x), for X G il. 

More precisely the following corollary holds 



(2.17) 
(2.18) 
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Corollary 2.4 Suppose that u <1 in O — \Jt^[r,T\^t, where r > and Q,t, for t e [r, T\, 
are smooth subdomains of Q and O is a smooth domain of x [r, T] and that u — 1 in 
Qt\0 then 

p{x,t) = —pc{u{x,t)) + constant{t), for all {x,t) e O 

and u satisfies 



Finally we remark that another form of the limit problem involves a parabolic equation, 
which is close to the standard Richards equation. Indeed if we set 4>{s) := Pc(0) — Pc{s) 
and denote by /3 the inverse function of 0, the function v :— (plu) is a weak solution of 
the problem 





in Qt, 

on on X (0,T), 
for X e n, 



with / e H{P{v)). 
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3 Existence of a solution of an approximate problem 
(S^) of Problem {S^') 

Let S be an arbitrary positive constant. In order to prove the existence of a solution of 
Problem (5''^) we introduce a sequence of regularized problems {Sg), namely 



(1 - U)t = dlv(^haiu)Vip + Pciu))^ + (l - rics)^Ss 

/ p{x, t)dx = 0, 

Vp.n = 0, 
"^{p + Pc{u))-n = 0, 
^ u{x,0) = Uo{x), 

where Uq, cg, ss and Sg are smooth functions such that Uq tends to uq in L^(Q) and cs, ss 
and Sg tend respectively to c, s and s in L'^{Qt), as 5 J, 0. In particular we suppose that 
there exists a positive constant C such that 



in n X (o,r), 


(3.1) 


in Q X (0,r), 


(3.2) 


for t e (0,T), 


(3.3) 


on dn X (0,T), 


(3.4) 


on on X (o,r), 


(3.5) 


for X e Jl, 


(3.6) 



sg>0, sg>0 and I si+ j s]<C. 

Jn Jn 

Moreover we suppose that Uq, cg satisfy 



< Urr, < Ua < 1 - S < 1 in fl 



and 



< Um < Cg < 1 — 5 < 1 in Qt- 
Adding up (3.1) and (3.2) we deduce the equation 

-div(^M^'{u)Vp + -ka{u)V{pc{u))^ ^sg-Sg- -f^{sg - Sg)dx. 



(3.7) 

(3.8) 
(3.9) 

(3.10) 



We formulate below an equivalent form of Problem (Sg). To that purpose we define the 
global pressure, V, by 



V -.^p+wiu) ^p+ r 

Jo 



kair) 



kair) + iikyjij) 



so that (3.10) gives 

-div(^Mi'{u)\/V^ ^ss-sg- -fjss - s^)dx. (3.11) 
We rewrite the equation (3.1) of Problem {Sg) as 

ut = Ar{u) + div(^r{u)M^^{u)VV) + r{cs)sg - nu)(^ss + -l^^(ss-ss)dx^, (3.12) 
where 

is continuous on [0, 1] and differentiable on [0, 1). Multiplying (3.11) by f^{ug) and adding 
the result to (3.12) we deduce that 



Ut = A^jj^'iu) + M''(u)Vr(u).VV + (^r(cs) - r(u)jss. 

This yields a problem equivalent to {Sg) , namely 

f Ut = Aip'^itj) + M''(u)V/''(^).VP + [ricg] - nu)]sg, in Q X (0, T), (3.14) 

-div(M^'{u)VV^ =sg- Sg- j-Jsg- sg)dx, inf]x(0,T), (3.15) 

^g^^ ^ J^V{x,t)dx ^ J^n^{u{x,t))dx, forie(0,r), (3.16) 

VV.n = 0, on9Qx(0,T), (3.17) 

V^^(M).n = 0, on dn x (0, T), (3.18) 

. M(a;,0) = M^(a;), iovxeQ. (3.19) 

In order to prove the existence of a smooth solution of {Sg), we introduce the set 

/C e C^+"'^(gr), Um<u<l-5}, 

where a e (0, 1), and we prove the following result 

Lemma 3.1 Assume {Hi) — {Hg) then there exists {ug,Vg) solution of {Sg) such that 

<eC2+"'^(Q^), Um<u''g<l-5 
andV^, VV^ e Ci+'*'^(gr) and AV^ e C°'t(gr). 
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Proof : Let be the map defined for all F e /C by T^{V) — W, where W is the unique 
solution of the elliptic problem, 



iQl 



V) 



-div[Mf'{V)VW^ ^ss-sg- -j^^{ss - ss)dx, in O x (0,r), 

J W{x,t)dx^ J n^'(V(x,t))dx, forte(0,r), 
vW.n = 0, on9fix(0,r). 

By standard theory of elliptic system (see [7] Theorem 3.2 p 137), we have that 

mZ^'"^ + iVW-lgt"''*' < dMqT"^ + D2. (3.20) 

For W solution of (Qy) we consider defined by T^{W) — V, where V is the solution 
of the parabolic problem, 

[ Vt = ^^^{V) + M^^{V)Vr{V).VW + [r{cs) - r{V)]ss, in Q x (0, T), 
iQlr) I Vij^{V).n = 0, on dn x (0, T), 

[ V{x,0)^uf^{x), inn. 

^Prom the standard theory of parabolic equations, we have that 



\V\q;' ^ <Ds[\W\q;' - +\VW\q:' ' )+D,. (3.21) 

Moreover defining by C the parabolic operator arising in {Q^) , namely 

/:{V){x,t) := Vt - ^^^{V) - M'^{V)VnV)VW - [/'^(c,) - nV)]ss, 

we remark that (2.8), the property (3.9) of cs and the fact that, by (3.7), ss is positive 
imply that 

jC{um) < and £(1 - 5) > 0. (3.22) 

Setting T :— T"^ o T^, the inequalities (3.22) ensure that T maps the convex set /C into 
itself. Moreover we deduce from (3.21) that T{}C) is relatively compact in /C. 
Next, wc check that T is continuous. Suppose that a sequence (Kn)mGAf converges to a 
limit V elCin C^+"'^(Qt), as m oo. Since (KOmew is bounded in C^^""' 2" (Qj,), it 
follows from (3.20) that the sequence {Wm '■= T^(Vrn))meN, where Wm is the solution of 
{Ql,J, is bounded in Ci+"'^(Q 2^), SO that as m ^ 00, Wm converges to the unique solu- 
tion W of Problem (Q^y) in C^+'^'^(Qt) for all (3 e (0, a). Moreover W G Ci+°'H^(gy). 
Purther it also follows from (3.20) that {VWm)m£N is bounded in C-*^"*""'^ (Qt), so 
that the solution Vn = T'^iWm) of Problem {QwJ is bounded in C^+'^^^iQ^). Since 
Vrn — T'^{Wm) — T{Vjn), {T{Vm))meN couvcrgcs to the unique solution V of Problem (Qw) 
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in C^+P'^iQj) for all (5 e (0, a), as m ^ oo, so that V = T'^{W) - o T\V). There- 
fore we have just proved that (T^ o T^{Vm))meN converges to o T^{V) in C2+^'^(Qr) 
for all (3 E (0, «), as m ^ cxd, which ensures the continuity of the map T. It follows from 
the Schauder fixed point theorem that there exists a solution {ug,Vg) of {Sg) such that 

< eKnC^+^'"-^{QT) and V^, V^', G C'+^''^{Qt), AV^ e C"'f (Qt). 



This concludes the proof of Lemma 3.1. Moreover we deduce from Lemma 3.1 the existence 
of a solution of (5"^), namely 

Corollary 3.2 Assume the hypotheses (Hi) — (Hg) then there exists {ug,pg) solution of 
(S^) such that < e C^+"'^(Qt), 

Um<u'^{x,t) <l-5 (3.23) 



4 A priori Estimates 

In view of (2.8) and (3.23) we deduce the following bounds 

= r(0)<r«(a;,t))<l = r(l), (4.1) 

= ^(0) <r(c,(x,t))<i = .r(i), (4.2) 

< K{u^) < k^iu'^six, t)) < kM = 1, (4.3) 

= ka{l) < ka{4{x, t)) < ka{Q) = 1, (4.4) 

< fejn„0 < Af^(Mn. (4.5) 

Pc(l) <Pc«(:r,t)) <p,(0), (4.6) 

Pc(l)-Pc(0) <7^^«(a;,t)) <0, (4.7) 

Pe(l)-Pe(0)<Q'^«(x,i))<0, (4.8) 

for all {x,t) e Q X (0,7). Next we state some essential a priori estimates. 

Lemma 4.1 Let {ug,pg) be a solution of Problem (Sg ). There exists a positive constant 
C , which only depends on Q, k^, ka and T such that 

r I ka{u^)\Vp^ + Vpc{u'^)\'dxdt < Cii, (4.9) 
JO Jo. 

r I \^I^5?dxdt < C, (4.10) 
Jo Jvi 
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and 



< - / / Vg{u'^).Vpc{u'^)dxdt < C, 



Jo Jn 



Jo Jn 

i-T 



Jo, 



\WCK)fdxdt < C, 
\Vg{u^s)fdxdt < C. 



(4.11) 
(4.12) 
(4.13) 



Proof : We first prove (4.9). Multiplying (3.11) by P = + 7^''(M^) and integrating the 
result on Qt = Q x (0, T) we obtain 

/ M'^«)|V(p^ + n^M)\' <\i (ss - ssf + h[ (p^ + 7^'^«))^ (4.14) 
for all h > 0. Moreover we have by Poincare-Wirtinger inequality that 

/ (p^ + 7^>^))^<cJ/ |v(p^ + 7^''«))r+(/ p^ + 7^'^«))'". 

JQt IJQt \JQt / J 

Using (3.3) and (4.7), it follows that 

/ (p^ + 7^^«))2 < c, [ |v(p^ + 7^-«))r + c,, 

JQt JQt 

k (u ) 

which we substitute into (4.14) with h = deduce, also in view of (3.7) and 
(4.5), that 

/ |V(p^ + 7^^«))r <C3and / |p^ + 7^^«)|2<C3. (4.15) 

JQt JQt 

Furthermore multiplying (3.1) by p^ and (3.2) by p^ -\-pc{u^).i adding up both results and 
integrating on Qt we obtain 

- / K)tPc«) + / ^.(OlVp^p + -fc.«)|Vp^ + Vpe«)r = /, (4.16) 

JQt JQt A* 

where 

We check below that first term on the left-hand-side of (4.16) and / are bounded. Denoting 
by Vc a primitive of Pc we have that 

/ p.K)K).=xri[n«)i. 
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Since Vc is continuous and it^ is bounded this gives 

/ PcK)K)tdxdt 
Moreover we have using (3.3) and (2.13) that 



< C4. 



- J^^ n^iu^s) [ncs)ss - nu's)ss + (1 - no) (/ - ss) 

+ J [(1 - ncs))ss - (1 - r «)) {ss + -f^iss - Ss))] [q'^M +Pc(0) 



dxdt 



(4.17) 



(4.18) 



dxdt. 



In view of (H^), (3.7), (4.1), (4.2), (4.7) and (4.8) we obtain 



This together with (4.15) yields / < C5C3 + Cq. Substituting this into (4.16) and also 
using (4.17) we obtain that 



Qt 



^-(OlVp^r + -A;a«)|Vp^ + Vp,{u^stdxdt < C7, 



(4.19) 



which implies (4.9). In view of (4.3), we also deduce from (4.19) the estimate (4.10). 
Next we prove (4.11). By the definition (2.9) of gi, we obtain from (3.10) that 

-div[M'^{u>'s)Vpi^ + ^A^«) = ss-ss- jj^ss - ss)dx. (4.20) 

Multiplying (4.20) by f^{ug) and subtracting the result from (3.1) we deduce that 

(<)t = ^/n<)A5(<)+c^^^(^«,(<)Vp^)-r(<)d^^(M'^(<)V(p'^))+^4r(c5)-/''(<)]. 



Moreover using the definition (2.7) of /'^ and Af^ we note that 



(4.21) 



dwfM'^«)/^«)Vp; 



div[k^,{u^)Vp', 
M^^{u^s)^{nu''s)).Vp>l + r«)(iz^;(M'^«)V(p'^) ), 



which we substitute into (4.21) to obtain 

{<)t ^ -rH)^g{u^5) - M^iu'^Mnu^)).vp'^ = uncs) - r«)]. (4.22) 
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We set 

D^{a):^p,{a)r{a)-Q^{a), (4.23) 

for all a e [0, 1], so that by the definition (2.11) of we have V^/^^) = Pc«)V(/''«)). 
Substituting this into (4.22), which we have multiplied by Pc{ug), we deduce that 

pc{u'sM)t - -^rK)Pc{<)^g{<) - M'^«)vD'^«).vp^ = PciOMncs) - nu^s)]- 

(4.24) 

Multiplying (4.20) by Df^{ug), adding the result to (4.24) and also using the fact that 



(i«t;(^iVF«)Z)''(w5)Vp^j = iVF«)VL"^«).Vp^ + Z}^«)(iit;(^iVF«)Vp^ ), 
we deduce that 



= Pc«)s5 



+ D^^iu'^s) [ss -ss- -f^iss - 55)) . (4.25) 
Integrating (4.25) on Qt and using the fact that the definition (4.23) of imphes 

Pc«)r(«^)-I^^(«^) = Q^«), 

we obtain 



where 



J := 



Pe«)«)tM - - / Q^{u^)Ag{u^)dxdt = J, 

/Jj J Qx 



I PcH)ss[ncs) - nu^)]dxdt 

JQt 

+ / (pe«)r «) - Qn<)) {ss -Ss- -f^iss - S,)^ dxdt. 



(4.26) 



It follows from (4.1), (4.2), (4.6), (4.8) and (3.7) that \ J\ < Cg. Substituting this into 
(4.26) and also using (4.17) we obtain that 



Furthermore we remark that 



< -- / VQ'^«).V(5«))cia;rft < Cg 

^■u)('^m) 



(4.27) 
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which together with (4.27) and the fact that VQ'*(?/^) = f''{u'^)Vpc{u'^) yields 

< - / Vpe«)V(^«))M < Cio. (4.28) 

J Qt 



By the definition (2.10) of C, we have -Vpc{u'^)Vg{u'^) = |VC«)p. This together with 
(4.28) imphes (4.11) and (4.12), which in view of (4.4) gives (4.13). This completes the 
proof of Lemma 4.1. 

In what follows we give estimates of differences of space translates of pg and g{ug). We 
set for r e R"^ sufficiently small: 

Qr = {xen, B{x,2r)cn}. 

Lemma 4.2 Let {ug,p'^) be a solution of Problem (5*^); there exists a positive constant 
C such that 



and 



JUr 

T 



p'i{x + ^,t)-p'^{x,t) 



{x,t)dxdt < C^^ 



(4.29) 
(4.30) 



/ / g{u'^){x + U)-9{0{x,t) dxdt<Ce, 

Jo JQr 

where ^eR^ and\^\< 2r. 

Proof : The inequalities (4.29) and (4.30) follow from (4.10) and (4.13) respectively. 

Next we estimate differences of time translates of g{u^)- 

Lemma 4.3 Let {ug,pg) be a solution of Problem {Sg) then there exists a positive con- 
stant C such that 

T"" / [g{u'^){x,t + r)-g{u'^){x,t)fdxdt < Cr, (4.31) 

Jo JQ, 

for all re (0,T). 
Proof : We set 

A{t) := / [giu''g){x,t + T)-g{u''g){x,t)]^dx. 



Since 51 is a non decreasing Lipschitz continuous function with the Lipschitz constant Cg 
we have that 



A{t) < Cg [g{u'^{x,t + T))-g{u'^{x,tmu'^{x,t + r)-u'^{x,t)]dx 



dx 



r rt+T 

< Cg / [giu^x, t + r))- g{u^{x, t))] / «)t(^, 0)dO 
Jn Ut 

L«(a:,t + T))-^«(x,t))l 
j\iJt I J 

div(k^(u'i)Vp^) + f^icsps - r «) [ss + f^iMy) - ss(y))dy 
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{x, 9)d9dx, 



where we have used (3.1). Integrating by parts this gives 

A(t) < C,y'^^ JJk^(u^)(x,e)VpUx,e)Vg(u^)(x,t + r) 

k^{u^s){x, e)Vp'^{x, e)Vg{u'^){x, t) dxde 



dxdd 



+ 



+ 



^«)(x,t + T) -g{u'^){x,t) 



K{x, t, T)dx 



(4.32) 



where 



K{x, t, t):= ^ [Ficsix, e))Mx) - r{us{x, 6)) \ss{x) + j^^iMv) - Si{y))dy 
Next we estimate the right hand side of (4.32). Using (4.3) we have that 



dO. 

(4.33) 



< 



fc^«)(x, ^) Vp^x, e)V^«)(x, t + r) dxde 

t JO, 

\ ( \T L ''"^^"^^ + \T L '^^^""^^ ^ + ^) i''^^^^) 

< li^l'^^ JjVp's(^,e)\'dxde + r JjVg(u^)(x,t + r)\'dxy 
Similarly we have that 

I k^{u''s)i^,e)Vp'^{x,e)Vg{u'^){x,t) dxde 
< ljVp^si^,e)\'dxde + T ljvg{u^){x,t)\'dxy 

Moreover using (4.1) and (4.2) we obtain from the definition (4.33) of K that 



(4.34) 



(4.35) 



ft+r r r "I r r 

\K{x,t,T)\< \ss\ + \ss\ + i- \s5-ss\dxde< \ss\ + \ss\ + f \ss-ss\dx 

Jt I J n JL J ci 



T. 



This together with (3.7) and the fact that the function g{ug) is bounded uniformly on /j, 
and 5 yields 



/ 9{us){x,t + T) - g{u'^){x,t) 
Jn L 



K{x,t,T] 



dx < Ct. 



(4.36) 



Substituting (4.34), (4.35) and (4.36) into (4.32) we deduce that 

A{t) < C,y^'^^ JjVp'^{x,e)\'dxde + yjVg{u^s)i^,t + T)\'dx 
+yjVg(u^s)M\'dx + Cry 
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which we integrate on [0, T — r] to obtain 

£ ^ A{t)dt < Cg(^J^ ^ J'^^ JjVp'^{x,9)\^dxd9dt + T JjVg{u'^)\^dxdt + CTT^ 

< Cgl^T JjVp'^{x,e)fdxde + T jjVg{u''s)fdxdt + CTTy 
In view of (4.10) and (4.13) we deduce (4.31), which completes the proof of Lemma 4.3. 

5 Convergence as 5 | 0. 

Letting 5 tend to 0, we deduce from the estimates given in Lemmas 4.1 and 4.2 the 
existence of a weak solution of Problem {S^). More precisely, we have the following 
result. 

Lemma 5.1 There exists a weak solution {u^,p'^) of Problem {S^), which satisfies 



r 

Jo Jn 



dxdt <Cn, (5.1) 

/' / \Wp^\^dxdt<C, (5.2) 
Jo Jci 

r I \Wg{u>')?dxdt < C, (5.3) 
Jo Jn 

r I [g{u^){x + i,t)-g{u^){x,t)fdxdt < Ci\ (5.4) 

where ^ G and \^\ < 2r. Moreover the following estimate of differences of time 
translates holds 

f [g{u''){x,t + T)- g{u''){x,t)fdxdt<CT, (5.5) 
Jo Jn 

for all re (0,T). 

Proof : We deduce from (4.10), (4.30) and (4.31) that there exist functions g^ and p^ and 
a subsequence ((ug^,p'^^))neN of weak solutions of Problem (Sg^) such that 

{9{usJ)neN tends to gf" strongly in L'^{Qt), (5.6) 
{PsJneN tends to weakly in L\0,T; H\Q)), 

as 5n tends to zero. Thus for a subsequence, which we denote again by 5„, we have that 

{g{ugJ)neN tends to g'^ for almost {x,t) G Qt- (5.7) 
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Using the fact that g is bijective we deduce that 

{ug )neN tends to := g~^{g^) strongly in L^{Qt) and almost everywhere in Qt, 

(5.8) 

as 5n tends to zero. Moreover we have in view of (4.13) and (5.6) that Vg{ugJ tends to 
'Vg{u^) weakly in LF'{Qt) as (5„ J, 0, so that by the definition (2.9) of g 

kaiugjVpciu'^J tends to ka{u^')Vpc{u^) weakly in L^{Qt) as 5„ i 0. (5.9) 

Letting Sn tend to in (3.23) we deduce that 

Um<U^{x,t) <1. (5.10) 

Moreover we deduce from (3.3) that 

/ p^{x,t)dx = 0, for almost every t e (0,T). (5-11) 
Multiplying (3.1) by </? e C, integrating by parts and letting 5„ tend to we obtain 

u^ifitdxdt ^ £ k^{u^)Vp^.V(fidxdt - £ (/^(c)s - f^{u^)sj(fidxdt 

— / uo(x)(p(x,0)dx, (5-12) 
Jn 

where we have used that Uq tends to Uq in L^(Jl) and that cs, ss and Sg tend respectively 
to c, s and s in L?'{Qt) as 5 J, 0. Similarly, multiplying (3.2) by </? e C, integrating by 
parts and letting 5n tend to we deduce that 



^ (l - u^^iptdxdt =-j^j^ [ka{u^)^jf + Vg{u^)^ .Vipdxdt 

~Io L (^^ ~ /''(c))^- (1 - nun)s)fdxdt - (l - Uo{x)y{x,0)dx, (5.13) 



which since 'S7g{u^^) — ka{u^)'Vpc{u^) coincides with (2.15). Next we prove (5.1). We first 
check that 

kaiu'^jVpl tends to kaiu'')Vp'' weakly in L^Qt), (5.14) 
as Sn tends to 0. Let ip e L'^{Qt), we have that 

Lr i^^^'"''-^^^- ~ ka{u^)Vp^y dxdt\ < + |/|J, (5.15) 

where 

II (kaiul) - k,{un)vpl<p dxdt 
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and 



II = j^^ K(u'^)^(vp^^ - Vp^^ dxdt. 



Using the fact that Vp^^ converges to Vp'^ weakly in L?{Qt) as 5n [ 0, we deduce, since 
ka{u^)^ e L'^iQr), that" 

|/|^| tends to as 5„ i 0. (5.16) 
Moreover we have by (4.10) that 



< 



< C 



(f K{ul)-kM) 

\ JQt 



^dxdt] 



1/2 



\Vpl\'dxdt] 



1/2 



(fi dxdtj 



X 1/2 



Since ka{ug^) — ka{u^) </?^ < 4(/7^ and since ka{ug^) tends to ka{u^) almost everywhere, 

we deduce from the Dominated Convergence Theorem that Ig^ tends to as (5„ J, 0. This 
with (5.16) imphes (5.14), which with (5.9) gives that 



tends to ka{u'^) 



weakly in L'^{Qt). (5.17) 



The functional v 



v^dxdt is convex and lower semi continuous from L'^{Qt) to R 



therefore it is also weakly l.s.c. (see [2] Corollary III. 8) and thus we deduce from (4.4), 
(4.9) and (5.17) that 



kaiun 



Vp^ + Vpciu") 



^dxdt < liminf^^^o / kMJ Vpl + ^PciuD 
JQt L J L 

< liminf^^^o/ ka{ul)\vpl + VpcK, 

JQt L 



dxdt 



dxdt 



which coincides with (5.1). Finally, we deduce respectively from (4.10), (4.13), (4.30) and 
(4.31) the estimates (5.2), (5.3), (5.4) and (5.5). This concludes the proof of Lemma 5.1. 



6 Convergence as | 0. 

The goal of this section is to prove Theorem 2.3. We first deduce from the estimates 
(5.2), (5.4) and (5.5) that there exists a couple of functions {u,p) and a subsequence 
{{u'^",p'^"))neN such that 

{u^")n£N tends to u strongly in L'^{Qt) and almost everywhere in Qt, 
lp^")neN tends to p weakly in 1.2(0, T; H\Q)), 
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as /In tends to zero. Moreover since 

there exists a function / G L^{Qt) with < / < 1 and a subsequence {f'^""' {u'^""''))nmeN 
of {f'^"{u'^"))n£N such that (/^"™ ))n„GAr tcuds to / wcakly in L'^{Qt) as tends 
to zero. Moreover we deduce respectively from (5.10) and (5.11) that < u < 1 and that 

J p{x, t)dx — 0, for almost every t e (0, T), 

which gives (2.18). As it is done in Section 6 in the proof of (5.1), one can first check that 

/ca(ii''"'")(Vp''"'" + Vpc{u'^""')) tends to ka{u){Vp + Vpc{u)) weakly in L'^{Qt), 

as i ^iid then deduce from (5.1) the estimate (2.17). Furthermore letting tends 
to zero into (5.12) we obtain, since lim^^^^o Z'^""' («) = x{^) for all s e [0, 1], that 

/ / uiftdxdt = / / kyi,{u)'Vp.'V(pdxdt — / / [x{c)s — fs](pdxdt 
Jo Jn Jo Jn Jo Jn \ / 

- [ Uo{xMx,0)dx, 
Jn 

which coincides with (2.16) and concludes the proof of Theorem 2.3. 

7 Numerical simulations 

7.1 The saturation equation and the numerical algorithm 

In this section we present numerical simulations in one space dimension. To that purpose 
we apply the finite volume method, which we present below. To begin with, we rewrite the 
equations (2.1) and (2.2) in the case that ft = (0, 1); this gives for {x,t) e (0, 1) x (0, T) 

< = d,(^k^{u^)d,p^^ + nc)s - riu^s, (7.1) 

(1 - u% = d,Q^ka{u'^)d,{p^+p,{u'^))^ + (1 - r{c))s - (1 - r{un)s. (7.2) 

Adding up both equations and using the boundary conditions (2.4) and (2.5) we obtain 

d^jf = - , , , ^ / , ^ dJpAu^)). (7.3) 



Substituting (7.3) into (7.1) yields 

^ + r{c)-s - r{u^)s. (7.4) 
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Ut = -dx 



' riun^^dxiPciun) 
1^ 



Moreover we deduce from (7.3) and the definition (2.12) of 7^'* that d^p^" = -9^(7^^(^i'')), 
so that in view of (2.3) we have 

p^ix.t) = -ni'iu^Mx.t) + n^iu^)iy,t)dy. (7.5) 

Jo 

In the sequel, we compare numerically the solution of (7.4) with the solution u of the 
limit equation in the case that u < 1, namely 

ut = -da: {kw{u)d^pc{u)^ + x(c)s. (7.6) 

We discretize the time evolution equation (7.4) together with the initial condition and 
the homogeneous Neumann boundary condition. The time explicit finite volume scheme 
is defined by the following equations in which /C > and J > Q denote respectively the 
time and the space step. 

(i) The discrete initial condition is given for i e {0, [1/^7]} by 

[U^t = u^{iJ,^)). (7.7) 

(ii) For i G {0, [l/J"]} and for n e {0, \T/]C\} the discrete equation is given by 



/c 

where 



(iii) For n e {0, [r/ZC]} the discrete Neumann condition is defined by 

[F% = and = 0. (7.9) 

The numerical scheme (7.7)-(7.9) allows to build an approximate solution, Uj^/c '■ [0, 1] x 
[0, T] ^ R for all i e {0, [1/J]} and aU n e {0, [T/IC]}, which is given by 

uj,k{x, t) = <, for all X e {ij, {i + 1)J] and for all t e (n/C, {n + 1)/C]. (7.10) 

In order to also compute the pressures, we propose the following discrete equation corre- 
sponding to (7.5) 

[pM]n ^ _7^([[/M]^) + j^flf^nip^"]]). (7.11) 

Finally, setting p'^{x, t) = p^{x, t) + Pc{u^^){x, t) we deduce that 

(ra)r = -7^([c/1^ +Pc([c/ir) + jM^^npr,). (7.12) 
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for all i G {0, [l/i/]} and all n G {0, [T//C]}. Similarly we propose a finite volume 
scheme corresponding to the equation (7.6), namely 



where 



n+1 



U:\=FlX,-F- + xiC:)S: 



i 1 



(7.13) 



for all (z, n) E {0, [1/J]} x {0, [T//C]}. 



7.2 Numerical tests 



For the numerical computation we take fi = 10~^, Pc{z) = 0, — z, ka{z) = (1 — z)^, 
kw{,z) = and s{z) = So{z), s{z) = 6i{z), where 6a is the Dirac function at the point a. 
Furthermore is given by the line with crosses, p'^ is given by the lines with diam and 
the limit function u corresponds to the continuous line. 

First test case: The case that c = 0, 7 and Mq = 1 on [0, 1]. We obtain at t = 0, 01 the 
following pictures 



0.0- ooooooooooooooooooooooooooooooooooooooo 



0.2 0.3 



o.a 0.9 



Figure 1 : t=0,01 

We note that, for /i small, the functions u and are very close. Here we only start 
with water and inject a mixture of water and air. The air immediately invades the whole 
domain. Figure 1 illustrates the result which we proved in this paper, namely that u'^ 
tends to the solution u of the limit equation (7.6) as fi tends to and moreover that the 
pressure is constant. This is indeed the case since m < 1. 



Second test case: The case that c = 0, 7 and uo{x) 



0, 1 on [0, 1/3] 



\ 0,7 on (1/3,1] ■ 
the following pictures for t = 0, 01 and for t = 0,1 respectively 



We obtain 
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Figure 2 : t=0,01 




Figure 3 : t=0,l 

The injection of a mixture of water and air (c = 0.7) takes place in a region of low water 
saturation. We first remark that both functions and u evolve very slowly. Here again 
we have that u{x,t) < 1 for all (x, t) e (0, 1) x (0, T) and we remark that the graphs of 
the two functions u'^ and u nearly coincide. 

Third test case: The case that c = 1 and uq(x) = < ^' i • We obtain the 

[ 0,7 on (1/3,1] 

following pictures for t = 0, 01 and for t = 0,1 respectively 
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Figure 4 : t=0,01 




<><><><>^^<><>0 0<><>0<> OOOOOOOO^^OOOOOOOO^^O 



0.0 0.1 0.2 0.3 0.4 0.5 O.S 0.7 O.B 0.9 1.0 

Figure 5 : t=0,l 

Here only water is injected; note that the saturation evolves rather fast. 
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